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ABSTRACT 



This report examines gas dynamic stabilization of high pressure, 
highly non-uniform electrical discharges for two electrode geometries. 

The geometries studied are the multiple-point to plane and the multiple- 
wire to plane. A successful analytical model is developed utilizing the 
anode extension concept, which describes the interactions and stabilizing 
effects of flow, separating the effects of convection and turbulence, in 
the point-plane geometry. Comparisons of theory and data are presented, 
illustrating actual observed and calculated breakdown voltages and 
currents for a number of cases. Flow stabilized high power density 
electric discharges are applicable to high power electric discharge lasers, 
magnetohydrodynamic power generation, and other new systems and technolo- 



gies. 



I . INTRODUCTION 



A. Background 

Electric discharges have been the subject of intensive study for many 
years. Numerous books, papers, reports, etc. have been published which 
characterize and analyse the various complex and interrelated processes 
which take place in discharges of several gas compositions, pressures, and 
electrode geometries. The primary objective in most of these studies has 
been to attempt to relate the atomic level processes to the observed 
macroscopic discharge properties such as breakdown voltage, current and 
overall discharge stability for use in engineering analysis and application 
in design. Some of the more established applications of this type of 
gaseous electronics include glow discharge lighting, the high voltage arc, 
corona discharge analysis around power lines and power line design, high 
voltage switching, spark plug and spark gap design, thermal printing 
devices and many others. 

More recent applications of this electric discharge technology have 
been generally aimed at systems which are geared to creating very high power 
densities within the discharge, thrusting physicists and engineers headlong 
into a direct confrontation with the stability problem, not just with the 
idea of “defining the envelope*', to use an aeronautical term, but to develop 
ways of expanding the natural discharge stability boundaries through any 
practical means. The recent high power applications of which we speak 
i ncl ude magnetohydrodynami c (MHD) power generation, plasma chemistry, 
electric discharge convection lasers. E-beam ion izer-sustainer lasers and 
others . 
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As a result of this growing requirement for increasingly higher power 
densities, many analyses have been done in the stability area and much 
insight has been gained in the many modes and regimes from which discharge 
instability can develop to cause a severe local breakdown in a delicate 
balance of parameters. One of the successful techniques which has appeared 
in the engineering of high power density systems is the pulsed system 
where-in the pulse time is shorter than the characteristic times of the 
various instability modes which must be avoided. The E-beam and pulsed 
ioni zer-sus tai ner lasers are examples of the success of this approach 
(Reference 1). However, present E-beam systems are large, expensive and 
delicate* In addition, in order to achieve high average power output 
from the pulsed systems, while the systems can be scaled up in size to any 
practical limit, once the size limit has been reached, the only other 
alternative is to increase the pulse repetition frequency (PRF) . But this 
option is also limited, usually by a characteristic optical cavity clearing 
time which depends on cycling out by flow the products of one pulse before 
initiating the next one. Failure to allow sufficient clearing time will 
result in cumulative space charge and gas heating instabilities with 
resultant low power breakdown. 

Due to the inherent limitations cited above and the implicit depen- 
dence of those systems on gasdynamics, the present research has been pur- 
sued in order to better define the role of turbulence in a high pressure 
discharge as has been observed by many researchers (References 2, 3, 5 

and 6). An understanding of the fundamental modes of gasdynamic coupling 
which have been observed to affect breakdown voltage and current may lead 



2 



to new ways to improve present high power density systems through bulk or 
electrode boundary layer stabilization techniques; and may also point the 
way for development of sufficiently high power self-sustained CW electric 
discharge systems which are inherently more compact, rugged, and less ex- 
pensive than the types mentioned above. 

B. Present Work 

Section A of Part II (Fundamentals of Gasdynamic Effects on a Highly 
Non-Uniform Electrical Discharge) of this report develops in detail a 
theoretical model for the calculation of breakdown voltage of a multi- 
point to plane (points positive) electric discharge based on the ‘*anode 
extension’’ concept initially proposed in Reference 7 and further discussed 
in Reference 3* The model incorporates flow and turbulence effects and is 
similar in physical basis to the “corona cloud” model of Hepworth, Klene, 
and Tozer discussed in Reference 8 for a static (non-flowing) non-uniform 
electric discharge. The present model successfully separates the gasdy- 
namics into the pure convective and pure turbulence effects, and isolates 
specific electrode geometry and gas parameters in individual terms. This 
facilitates understanding of the fundamental relationships between the gas 
properties, electrode orientation and geometry, and the individual gasdy- 
namic coupling mechanisms. 

This section also presents a qualitative discussion of the breakdown 
characteristics of the uniform and non-uniform discharge geometries high- 
lighting the similarities in a qualitative and semi -quant i tat i ve manner. 

It is shown that Meek’s breakdown criterion for the non-flowing, uniform 
field and the equation developed in Section A for the non-uniform, flowing 
discharge transform into very nearly Identical expressions in the middle 
ground of a flowing uniform field. This illustrates the general nature of 
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A 

the equations of Section A and provides straightforward physical insight 
into the final breakdown mechanisms in the geometric extremes. Radial 
diffusion and spreading is considered as related to the mean or "drift'* 
field in a highly non- tni form geometry. 

Section B of the report discusses the typical discharge section used 
in Section A with the three rows of pins replaced by three thin wires (see 
Fig. lb). A preliminary analysis which was done in the terms of Reference 
9 basically indicated the potential for improving discharge input power by 
a factor of 2-3 with the wire anode. Various degrees of discharge input 
power enhancement have been noted in similar configuration modifications 
in connection with MHD laser/power generation applications (Reference 10). 

Recent experimental results in this area are presented and discussed 
in relation to the points-plane configuration. To date due to power 
supply limitations no improvement in performance has been noted, and some 
recommendations for further research in this area are presented, when a 
larger power supply is available. 

Section C of Part I! discusses electrode design geometry and electric 
field distributions for optimum electron production and for best reduced 
electric field strength (E/N or E/P) for laser pumping. 

Part III presents a data summary of current and past point-plane and 
wire-plane laboratory work, reviews the theory and experimental trends, 
indicating important correlations and features of the anode extension model. 
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II. FUNDAMENTALS OF GAS DYNAMIC EFFECTS ON A HIGHLY 
NON-UNIFORM ELECTRICAL DISCHARGE 



A. Prediction of Flow and Turbulence Effects on Breakdown Voltage 

I . Description of Discharge Geometry for Analytical Formulation 

A significant body of breakdown data has been generated in var- 
ious highly non-uniform electric field geometries with positive and neg- 
ative polarities at the highly stressed electrode. The poi nt- to-pl ane 
geometry with the points positive, however, has been the particular case 
through which the present analytical model has been developed; but the 
model is sufficiently general that the basic flow and electric field 
interactions described should also be useful in analysis, interpreta- 
tion, and trouble shooting of many other non-uniform geometries and 
systems both with and without flow. 

The specific apparatus used in the current research, described 
in detail in Reference II consists of multi-point to plane electrode 
arrangement (see Fig. la) which is mounted in a wind tunnel with avail- 
able test section velocities on the order of 200 m/s of ambient pressure 
air. Varying degrees of turbulence can be introduced by placing turbu- 
lence generating screens upstream of the electrode set-up. Two flow 
modes are discussed in this report, namely, 

a. Normal Flow - directed through the turbulence screen, past 
the positive multi-pin anode, through the interelectrode gap, and thence 
through the grounded cathode, which is designed to approximate a plane 
with low blockage to permit maximum gas flow rates (see Fig. 2a). 

b. Reverse Flow - same as above except that the anode and cat- 
hode positions are reversed (see Fig. 2b) and the turbulence screen was 
omi tted. 
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While our earlier effort which included cross-flow tests (Fig. 
2c) was based on an assortment of turbulence screens which produced a 
variety of spectra and intensities, the present work (which includes 
verification of past trends, and additional data taking for better curve 

definition at low speeds) is done utilizing only screen #3 design. This 

screen proved to be most successful in the prior work (References 3, 7, 

ID- 

G^rieral Description of the Discharge 

With no voltage applied across the gap, and no gas flow, there 

exist a small number of random electrons present in the gap which are 

produced by outside sources (Reference 12). As voltage is applied these 
electrons under the influence of the resultant electric field begin to 
move toward the anode producing a small saturation current whose magni- 
tude is limited by the number of electrons present and their average 
drift velocity in the field. In a highly asymmetric geometry such as 
the point-to-plane, the electric field produced is very non-uniform, 
being highest at the point surface and falling off rapidly in a hyper- 
bolic fashion with distance from the point as illustrated by the axial 
field distribution in a hyperboloi dal point-to-plane gap (Reference 17) 
which is assumed in the following analysis: 

R = point radi us 

x= axial distance from point 

Ep = field strength at point = Riog ^ td/k) 



6 



V ^ app lied vol tage 
d 5 gap length 

Due to this field asymmetry, as the voltage increases further, the field 
strength at the point soon reaches a value which permits the free elec- 
trons moving toward the point to gain enough energy to ionize by colli- 
sion some of the gas atoms and molecules near the point surface, thus 
additional electrons and positive ions are produced resulting in a rise 
in the current. The voltage at which these first few ionization events 
occur is called the corona initiation voltage and is usually designated 
, i.e., when £p - ^ ^ • 

2V 

Rlog(4d/R) 

E^Rlog(4d/R) 



(2j 

(3) 



From Reference 13 the threshold field can be calculated from Peek's 

Law ( 1929 ) modified by Gary et al (1972) for a point-to-plane type 
geometry, i .e. , 



E = E MA 
c o 



1 + 




where 



(4) 



E = 31 KV/cm (near spark breakdown of air) 
o 

M = coefficient describing conductor surface (0. 6-1.0) 

K = 0.308 cm^ 

R = radi us of point 

A 5 0.392 P/(273+T) = 1 for air at STP 
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P = pressure in torr 
T = temperature i n C 

Continued increases in voltage moves the high field region further and 
further into the gap. At the position where the reduced field value 
(E/P) is approximately equal to that required for ionization by elec- 
tron collision (E./P) , free electrons will begin to produce avalanches 
into the positive point. From this critical distance , each 

avalanche will produce an exponential increase in electrons and positive 
ions left in the aftermath. The number of electrons produced in the non- 
uniform field is given by the well known expression (Reference 12) 

Xz 

/ adx 
^1 

n = n^e (5) 

where 

n = original number of electrons 
o ^ 

a = Townsend's first ionization coefficient 
X = avalanche travel distance 

The coefficient a is a strong function of the electric field as illus- 
trated by the usual semi -analyti cal expression found in many references - 

BP 

a = (6) 

where A and B are gas constants, P is pressure. It is necessary, 
therefore, to integrate the exponent since in the region under 
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discussion E and therefore a vary greatly with x as the avalanche 
rushes toward the point, Fig. 3* An average or effective ionization 
coefficient can be defined as follows: 



X 

a = 0 . dx (7) 

o 

X 

o 



which will produce the same number of particles in traveling the same 



distance, and then we can write 



a(x-x ) 
o 

n = n e 
o 

similar to the simpler uniform field case. In the following discussion 
the bar wil 1 be dropped for convenience, and x^ assuned to be 0. Due to 
the much greater electron drift velocities, the electrons from these 
avalanches are quickly drawn into the positive point leaving the positive 
ions behind in a maelstrom of activity including electrostatic repulsion 
from the positive point, diffusion in ail directions, neutralization by 
recombination with other electrons or charge transfer from negative ions, 
etc. All of these processes tend to dissipate the space charge concen- 
trations and smooth the non-uniformity of the electric field which is 
generated by the dense space charge near the point, and by the point 
i tsel f . 

However, if the applied voltage is raised high enough (near 
breakdown), these dissipative processes lag behind production processes, 
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and the space charge density and field intensity can become of the same 
order of magnitude as the applied field. When this happens, the portion 
of the gap external to the high field region is no longer able to dis- 
tinguish the applied field nor the position of the point, and the dense 
space charge actually assumes the role of the point forming a "virtual" 
anode (Reference \k) effectively diminishing the gap as seen from the 
cathode. The effect of this "extension" of the point is to require 
dissipation of the applied voltage over a shorter gap, thus raising the 
"drift" (low field) portion of the gap to a uniformly higher level (Fig. 

4). Reference 1^ illustrates how, in a non-uniform geometry, the presence 
of space charge can modify an originally hyperbolic field distribution 
to produce a uniform field over most of the gap in concentric cylinders 
where the inner electrode radius is much less than the outer. 

Another important effect which must be accounted for is the 
radial diffusion of the space charge in the ionization region which 
essentially distorts the shape of the original point, making the resultant 
virtual anode appear extended in the radial as well as the axial dimen- 
sion (see Fig. 5) . 

It will be seen in the following analysis, that accounting for 
the two effects (shortened gap, and modified radius) provides much in- 
sight into various discharge phenomena heretofore not well explained 
and leads to a very general form of the breakdown equation for a highly 
non-uniform geometry in a turbulent flow field. The breakdown equation 
is derived by accounting for modifications by the flow of the quasi- 
static space charge distributions which develop in the non-flowing run- 
ning discharge case« The trends revealed and results obtained are 
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fundamental and adaptable to different gas mixtures, geometries, and flow 
field parameters,, 

3 * Analytical Development of the Breakdown Equation with Flow 
Interactions 

The exponential charged particle multiplication process in high 
reduced fields has already been mentioned. Since the current is produced 
by the motion of these particles between the electrodes, the current 
density can be written as a summation of the individual species contribu- 
tions. For the present it is assumed that the positive ion and electron 
currents dominate the discharge (Reference 8). The current density will 
then be given by 



j 



en . V . + en V 
i i e e 



( 8 ) 



where 



j = local current density 

e = electronic charge (1.6x10 coulombs) 

V. = ion/electron drift velocities 
i ,e 

n. ^ E local charge densities 



Considering only the high field ionization region, if each elec- 
tron entering causes an electron avalanche, the average number of elec- 
trons produced is equal to the multiplication factor plus one, while the 
number of positive ions produced would be just the multiplication factor. 
The multiplication factor is usually so high, that the number of positive 
ions and electrons created in this region are identical. Therefore 
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( 9 ) 



n. - n - n 
1 e 

Due to the relatively low mobility of the positive ions, the 
majority of electrons continue into the anode and a region of growing 
positive space charge is produced in the immediate vicinity which will 
slowly begin to drift away due to electrostatic repulsion from the anode. 
Under these conditions the current density can be estimated by the 
following 



j = en (vj+v^) (10) 

Under conditions of low positive space charge, v^ is much greater than 
V. and the equation reduces to 

j - env^ 

But if the positive space charge field becomes of the same order of mag- 
nitude as the local applied field, electrons may be retarded or even 
trapped by the space charge cloud. In this case the electron drift 
velocities may be reduced to the same order of magnitude as the ion 
velocities, or to zero and even negative if the electron is trapped and 
dragged in the opposite direction by high positive ion fields. In this 
case, the current density approaches 

J - envj (11) 

and the particle density may then vary between 
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n - 



to n - 



( 11 ) 



ev ev. 

e I 

It is now proper to return to the "anode extension" concept introduced 
above. As indicated, when the space charge fieid near the anode 
approaches the same order of magnitude as the applied field, flux lines 
from charges within the remainder of the gap and from the cathode will 
no longer be able to reach the actual point, but will instead terminate 
on charges in front of the point. This essentially shortens the gap and 
to a first approximation, results in dissipation of the applied potential 
over a reduced distance, raising the mean electric field over the 
remainder of the gap (Fig* 5). 

In order to continue the development, an estimate of the length 
of the anode extension region (hereafter designated by a (S ) i s 
required. A reasonable estimate of the beginning of the region of dense 
positive charge would be the point at which the applied electric field 
approaches the breakdown value for the type of gas mixture involved. In 
our example of atmospheric air, it is assumed that the anode extension 
in a non-flowing discharge begins at the value E(x) = 25,000 volts/cm 
(E/P - 35). This assures that the majority of electrons approaching 
this position avalanche toward the positive point, depositing the slower 
moving positive ions to form the anode extension (see Fig. 5). 

As soon as a space charge concentration develops, however, there 
are numerous forces and processes which act on and within it to dissipate 
it. Among the important process to be considered are: 

(a) diffusion in all directions. 



13 



(b) electrostatic repulsion by the electrode of like sign, in 
this case the nearby point, 

(c) attraction by the opposite electrode which may be 
insignificant in this region in very non-uniform fields, 

(d) recombination with supplementary electrons supplied from 
various sources, i.e<,, photoelectrons, cathode emission, 
free electrons, detachment from negative ions, etc. 

In a "running*' discharge, i.e,, one which has reached a pseudo-steady 
mode of operation, a quasi-equilibrium condition is ultimately reached 
in which the charge production and dissipation mechanisms locally balance 
each other. The result of this balance is a quasi -steady anode extension 
with a finite and calculable extent; a distorted radial dimension of the 
virtual anode due to radial diffusion, and an increased mean electric 
field in the drift region of the discharge (see Fig, 5). Reference 3 
alluded to these two important dimensional changes. 

To calculate the equilibrium boundary of 6 in the axial 
dimension between anode and cathode, electrostatic repulsion, convection 
(in the flowing case) and recombination must be balanced, diffusion being 
small in the high field region compared to the axial component of ion 
drift. The radial boundary is determined by radial diffusion, and re- 
combination, and in the case of crowded multi-pin electrodes, by symmetry 
consi derations, 

a. The Anode Extension - Axial Dimension - 6 

Since the length of 6 is determined primarily by a balance 
between the ion drift out of the dense space charge region, and the dis- 
sipation of that charge by recombination, an expression for recombination 
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is required. The form assumed here as in Reference 3 is the three-body 
solution of Hinnov and Hirschberg, namely that the recombination time 



t 



REC 




( 12 ) 



where 



t 



p = recombination coefficient 

n = particle density 

« = characteristic recombination time 



Now, the positive space charge created in the vicinity of the 
anode will be subject to electrostatic repulsion until it recombines with 
an electron. Therefore, the travel distance of an individual ion during 
its 1 i fet ime wi 1 1 be 



X 



d 



■REC 



(13) 



where v. = k.E = ion drift velocity. Inside this drift distance will 

I I 

be much positive space charge, and outside it, most of the charge will 
have been neutralized by recombination. Therefore, since 



X . 2: 6 
Cl 



comb ining Equations 12, 13, and 14 



« = V. (^) 

2pn 



Substituting for n from Equation 11 



V. ev 2 V. ev. 2 



(u) 



(15) 



( 16 ) 
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So far, flow has not been introduced so the above equation is 
general. Taking flow into account, the following notation will be used: 



5^ = anode extension in flowing gas 

5 = anode extension with no flow 

o 



5 

5* = flow value normalized by no flow value (^) 

o 

4 

Typical flow velocities in practice range up to about 1-2x10 cm/ sec. 
Typical ion drift velocities in atmospheric air at the outer limits of 
the anode extension (E/P = 25,000/760 - 33) also run in the neighborhood 
of 5“8x10 cm/sec. 

Defining the slip parameter S = u/v. = , significant flow 

I K . tl 
I 

effects can be expected when the slip parameter is of order one, especial 
ly in a dense collision dominated discharge where the mean neutral gas 
motion is easily transfered to the ion population. As can be seen in the 
sped f i c examp le , 



S 



u ^ 1x10^ 2x10^ ^ 

k.E 8x10^ 5x10^ 




(17) 



and substantial space charge redistribution should be anticipated. 

Going back to Equation 16, and rewriting to include flow effects 

we have 

a. No electron/ion coupling: 
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and 



5 

u 



(v.+u) 




■V 

e 



2 ( 1 _) 



k.E(l+S) 

_j 

. 2 



2 




(S 

o 



V. V 

I e 

.2 

J 




(k.E)v ^ 
I e 



,2 

J 





(1+S) 



( 18 ) 



(18a) 



(18b) 



b. Heavy electron/ ion coupling (high space charge fields): 



5 

u 



Ju Ju 



and 



where 



(S 

o 



V 



3 



J 






(k.E) 

.2 

J 




X _ (HS)^ 
J '■ 




(19) 



(19a) 



(19b) 



From these considerations it can be seen 
velocity term (l+S) will vary with space charge 
and we will write a generalization of l8b and 19b 



that the power of the 
density (or current) 
as fol lows : 



5* = 0+S)'" 



( 20 ) 
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where m Is a variable exponent. 

Since continuity considerations require that the total current 
passing any equi potent ial surface within the discharge be constant the 
total current can be represented by the electron current entering the 
anode where ion contributions are zero, so 



j 




. af(6)5 

J ^ 

■^o 



( 21 ) 



where is the electron current flowing into the ionization region 

(defined by which is proportional to 6 ) , the exponential is 

the multiplication factor, and f is an as yet unspecified function of 
6 , therefore - 

a f5 

with flow: j = j e ^ (22a) 

a f5 

no flow: j = j e ° ° (22b) 

•^o •'oo 

fa 6 (a*6*-l) 

then j* = j^*e ° ° (23) 



NOTE: The * value of a quantity is used consistently to indicate the 

value of that quantity in a flow field normalized by its no 

flow value (except for pressure - P* to be introduced later is 

equal to P/P = local value/value at STP = 760 TORR) . 
o 
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The function f(6) has been introduced to relate the critical 

ionization distance (that point at E(x) = E^ppL , ED ^SPACE CHARGE 
Ej) to the anode extension (S « It is apparent that 6 is proportional 
to i«e«, as 6 increases either due to convection or increasing 

applied voltage, also increases due to increased applied or dis- 

placed space charge field, and conversely, as the anode extension shrinks, 
will decreaseo However, due to the non-linearity of the fields 
involved, a variation in 5 will not produce a linear change in x.^ 

To account for the non-linearity, f(6) is assumed to obey a simple 
power law form 

f = C6*" (24) 

Squaring j* and substituting into Equation 20, yields the following: 



j^“ exp[(2C6""c«^«^(o<-6*-l)] 

Rearranging and taking the natural log of both sides, 

£n6* + 2 Ca 6 5 ^(a*6*-l) = miln(l+S) - 2nj ^ 
o o o 

and solving for fi" and collecting (S* terms on the left side, 



(25) 



*n-l 

X ^ * 



' .):^ Jln5* = ( 



2cct 6 a 
o o 



( 1~^S) - ■2'nj 

ZCa c a^' o 



( 26 ) 



o o 
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This is the basic anode extension equation in a flow field, where again 
the * values are the flow values compared to the no flow values, i.e,: 



5 . = 5 ^/ 5 ^ 



a* = a /a 
u o 



S = u/v. 



Jo* = ■icu''Joo 



= normalized anode extension 
= normalized ionization coefficient 
= the ion slip parameter 
= normalized saturation current 



C = CONSTANT 



The exponent n of 5* turns out to be a free parameter of 
the problem whose exact nature is not yet well established, but is under 
investigation. The parameter m is the power of the velocity term 
(1+S) which analytically was determined to range from approximately 1 
to 3 , although experimental data indicate that occasionally a higher 

value is required. The reasons for this disparity are also under active 
Investigation. 

Inspite of the limitations cited above, this equation, even in 
its preliminary form (Equation 26) describes all the previously observed 
breakdown voltage characteristics on the particular apparatus used in 
the initial tests for various velocities, turbulence levels, and flow ori- 
entations tested, even including some seemingly anomalous points in 
earlier tests which radically departed from the norm as reported in 
References 3 and 7^ The equation also predicts some salient features in 
the low velocity ranges which had not previously been observed due to 
strong emphasis on high speed testing. These low velocity features have 
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been confirmed in recent tests and will be described in detail in part 6 
of this section of the report. 

To check the mathematical correctness of the form of this equa- 
tion, when there is no flow, S = 0 , and by definition, the (*) terms = 

1 . Therefore, the three In terms go to zero, a* = 1 and we are 

left with the identity 0 = 0 as i t should. 

The peculiar geometric factors of a problem are all contained in 

the 6 terms. Gas and electric field parameters enter through the a 

and S terms; and turbulence effects enter through a* as will be 

explained in the following sections. Secondary electron production 

processes are included in j * . 

o 

To obtain a solution to Equation 26, requires some auxilliary 
equations which are developed below. 

Auxilliary Equations 

a. The no-flow anode extension (6 ): 

o 

The basic assumption of the anode extension model is that 
the gap breakdown potential is a linear function of the effective gap 
length, i.e., the actual gap length less the anode extension (which in 
effect is a portion of the gap already in a state of local breakdown). 
There fore: 



Vb ::: d" :: d - 5 = E(d-6^) + (27) 

o o o c 

Vb = no flow breakdown voltage 

d"* = effective gap 
d = actual gap length 
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5 = no flow anode extension 

o 

E - mean electric field in effective gap 
= corona initiation voltage 



For a specific geometry, Vb^ can be calculated from the condition that 
E(x) = Ej at X = , resulting in 




E. D E 

~ — L - ^ ^ = 

P ~ (R+x) P 



D 1 2Vb 

R 1 o 

(R+C.6^j "F Rlog(4d/R) 



(28) 



and 



E. 

Vbo = ■^P*(R+C-5^)log(4d/R) (29) 



for the hyperboloi dal point-to-plane. 

An analytical form can also be determined for as in 

Section 2 by assuming that is a linear function of the point field, 

i,e., when the region in the vicinity of the point reaches the reduced 
field strength for ionization. Again for the hyperboloi dal point, this 
results in 



E 2V ^ E. 

= PRlog( WR)" ~ ^ 

and 

C E. 

s -^P*Rlog(4d/R) 

but Vb - V = E(d-5 ) 

o c o 



(30) 



(31) 
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and inserting Equations 29 and 31 into 27 yields: 



5 c 
o 



log(l»d/R) + 1 



( 32 ) 



b. Ionization Coefficient Equati 



(l) No Flow - use Equation 6 with 0 subscript 



a = A Fe 
o 



‘e 



(33) 



(2) With Flow/Turbulence Effects 



The primary coupling of a with the flow is result of 
pressure/density fluctuations in a highly turbulent flow. Therefore, by 
substituting P ± AP into the above equation and taking a simple average 
[a(P+AP) + a(P-AP)j /2 , and eliminating second order terms results in - 



a - 
u 



or 



a 



a* = 



a 

o 


BAP 


E 


2 


e 












’ BAP 


1 


■ E 


2 


e 



+ e 



BAP 

E 



( 3 ^) 



BAP 



+ e 



(35) 



The pressure fluctuations due to turbulence can be related to the mean 
flow velocity through the local pressure-density relation in an incom- 
pressible flow, i.e,, the total pressure equals the sum of the static 
and dynamic pressures. 
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( 36 ) 



1 2 

P = Po ^ ^ pu 



Assuming local test section total pressure remains constant in the 
presence of turbulence, and that the fluctuating component of 

velocity, is proportional to U , the mean flow velocity, then 



P - P = AP 2 5pu 
o 



.1 



(37) 



The quanti ty ( VlT^/U) is called the turbulence intensity (Reference 
18). At constant temperature, p ~ P ; therefore for constant T , 



AP ^ CpU^ 



(38) 



but also 



S = U/k.E U = Sk.E 



u/k.E 



u = U = Sk.E 



k. = 



o P 



and 



AP = 



?P(Sk,E)' 



(39) 



Substituting into Equation 35 yields finally - 



a« = 



S^E 

e +e 



2 2 

B?P k S E 
o o 



m 



which when u = 0 = S = 0 , reduces to a* = 1 as it should. 
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c. Breakdown Voltage Equation with Flow and Turbulence 

The ultimate objective of the above derivations is to arrive 
at an expression for computing the fundamental gas dynamic effects on the 
breakdown voltage of an electric discharge. Approximating the breakdown 
voltage ratio as 



Vb d-6 

Vb* = from Equation 27 , (4l) 

o o 

then in final form, 



5 

1 - (^) 5 * 

Vb* ■ ; (W) 

(1 ->■) 

which also reduces to Vb* = 1 when U = 0 6* = 1 . 

It should also be noted that when 6* 0, Vb*) , the 

max 6 ’ 

same relation which was postulated in Reference 3 based on the anode 
ex ten s i on con cep t . 

d. Summary of Equations 

For convenience of reference in the following discussion, 
the basic breakdown potential and auxilliary equations are repeated 
be 1 ow : 
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5 

(l) Normalized Anode Extension 5* = ^ : 



5' 



% 






(26) 



(2) No Flow Anode Extension 5 : 

o 



P*(|-|^)log(4d/R) + 1 



(32) 



(3) No Flow Ionization Coefficient 

BP 

a = 
o 



(33) 



a 

(4) Normalized Ionization Coefficient cx*= — : 

ct_ 



2 2 

B?P k S E 
o o 



2 2 

-B?P k S E 
o o 



a* = TT 



+e 



(^ 0 ) 



Vb 

(5) Normalized Breakdown Voltage Vb* = 

o 



Vb* 




(41) 
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5. Qualitative Discussion of the Equations 

The key to a clear understanding of the equations lies in a 
simple analysis of the normalized anode extension 6* , The basic form 
of this equation is 






In fact it was derived from Equation 20. 



5* = 



_ (HS)' 



..2 



2 2 

But S ^ u and j* n* which gives 



m 

n* 



m 



m 



which is essentially the relationship derived in Reference 3. Thus, the 
anode extension is a product of two factors, the pure convective velocity 
which, in normal flow, tends to lengthen 5* , and the inverse of the 
charged particle density, which tends to shorten 6* due to faster 
recombination at higher n* and consequently results in less particle 
dri ft. 

In the reverse flow mode, the 6* equation can be approximated 

by: 



6 * 



(l+s)'" 



(45) 
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where u and therefore S are negative, which decreases 6* in a rapid 
fashion with increasing negative velocity. It is interesting to note that 
in the reverse flow approximation, a* cancels out and Vb* is not 
influenced by turbulence. 

The breakdown voltage behavior is opposite in sense to 5* , so 
when 5* increases, Vb* decreases and vice-versa. Note: frequent 
reference to Figure 6 will be helpful in interpreting the following. As 
mentioned in the reverse flow mode, 6* is quickly reduced and therefore, 
Vb* immediately rises (Part 2 of Figure 6) and quickly approaches an 
asymptotic limit (Part 1 of Figure 6) as 6* 0 corresponding to the 

requirement to breakdown the entire physical gap. 

In the normal flow mode (flow from anode to cathode), the initial 
low gas velocity tends to lengthen 6* , which causes a decrease in break- 
down voltage (Part k of Figure 6), but due to the exponential dependence 
of n* , j* , and a* on U (through S) , with effective turbulence the 
a* terms in the denominator of Equation 26 soon prevail, causing 6* to 
decrease and Vb* to increase at higher velocities and turbulence levels 
(Part 5 of Figure 6). The velocity at which this transition occurs is a 
direct function of the effectiveness of the turbulence introduced. If the 
turbulence is inhomogeneous or ineffective, the Vb* curve may not reach 
a minimum in which case, further increases in velocity may actually reduce 
power input to the discharge (Part 7 of Figure 6, dashed curve extension). 

It is interesting to note also that in the normal flow mode with 
effective turbulence, after initially increasing, 5* can again approach 
0 as its absolute limit (Part 6 of Figure 6) due to the high recombination 
rate, thereby confining the normal flow Vb* rise to the same asymtotic 
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limit as applies in the reverse flov/ case, and again cor respond i ng to the 
requirement to breakdown nearly the entire physical gap. 

It should also be noted that this asymptotic behavior, which has 
been observed in many cases, (References 2 and 3), poses a fundamental 
limitation on the stabilization of the discharge potential by gas dynamic 
means. The final breakdown may possibly be the result of either of two 
separate sequences: 

(a) the relatively uniform external “drift" field may reach 
such a high value that a random electron avalanche originating in the 
vicinity of the cathode can take enough energy from the field to com- 
pletely cross the i nterel ectrode gap, resulting in an internal break- 
down of the gas, or, 

(b) if the cathode material has a low work function it may 
begin emitting numerous secondary electrons prior to the drift field 
reaching the gas breakdown potential in which case the breakdown proba- 
bility will rise and a cathode limited breakdown may occur. 

As turbulence intensities are increased for a given gap length, 
the reverse flow region is virtually unaffected, except to hasten the 
rise of Vb* to the asymptote, while in the normal flow direction, 
increasing turbulence results in higher a* values with an eventual 
upswing in Vb* due to the resultant 6* reduction. Figure 7 illus- 
trates the progressive modification of the Vb* curves with increasing 
turbulence. The asymptotic behavior observed at high normal velocities 
and turbulence levels confirms the apparent rise in the Vb* curve and 
explains logically the apparently anomalous “short gap, high current" 
data points reported in References 3 and 7. 
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The analytic solutions illustrated in Figure 7 also provided an 
additional area for test and verification of the breakdown equations. It 
can be seen from the Figure that even in the "short gap, high current" 
cases, there should be at least a small region of decreasing Vb* at very 
low normal flow velocities due to the non- turbulent convective effect. 
Recent low speed tests have confirmed this characteristic as the data in 
Figure 8 show. 

Another qualitative assessment of the equations can be made by 
observing the behavior of the curves as a function of the gap length. 

Figure 9 illustrates this behavior for a 1 and 5 centimeter gap which 
encompasses the ranges tested. The trends of the curve shift in both 
normal and reverse flow with changing gap length were as indicated by the 
theoret i cal curves. 

Section 6 provides a more quantitative comparison of observed 
data and the appropriate theoretical curves. However, before moving into 
that subject, a very interesting and rather striking qualitative comparison 
can be made of the form of Equation 26 and Meek's Breakdown Criterion 
(as presented in Reference 12) for the uniform field in a static (non- 
flowing) discharge. As will be seen, despite the very different conditions 
in the two discharge situations, and the very different approaches used in 
deriving the breakdown criteria, the resulting breakdown equations are 
remarkably similar. 

Since a description of the non-uniform field anode extension has 
already been developed, a short description of the Meek-breakdown scenario 
in the uniform field will help put the remainder of this section in 
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perspective. The Meek Criterion is based on the relative space charge 
field strength at the head of an electron avalanche compared to the applied 
electric field. If the avalanche travels a certain critical distance under 
the influence of the applied field, the space charge field at the surface 
of the space charge region will become the same order of magnitude as the 
applied field. When the field ratio K = “^0,1-1, corresponding to the 

avalanche traveling the critical distance , the space charge field can 

no longer be ignored. The Meek Criterion from Reference 12 is written as 
fol lows : 

ax + £n— = 14.46 + £n*!^^ -O.Silnpx + iln x (46) 

P pa 

The method of solution and some results are discussed in detail 
in the reference. What we will concentrate on are some straight forward 
substitutions into the above equation for comparison with the anode 
extension model. If the gap length is adjusted such that an avalanche 
becomes critical Just as it reaches the anode, then the electric field 
lines will be distorted in the direction of the space charge, as in 
Figure 10a and the concentration of charge will produce a non-uniform 
field and, in fact, just as the avalanche touches the anode, it acts as a 
small area (like a point) of anode extension. The axial length of this 
extension is 2r , but r (the avalanche radius) is proportional to the 
distance it has traveled, so from 2r 6 x , 

X ~ c5 (47) 
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and it can already be seen that the situation looks very much like the 
point-to-pl ane geometry. After making this substitution, and then 
rearranging and substituting non-dimensional i zed flow variables for all 
the quantitites, the result is 



6" - (— f— )— ^n6" = ~ ( V — 

ca 6 a* ca 5 a 

o o o o 



11 o 1 “ a*kp-:^ 

)^nS + (— L^)i_ - 1) (48) 

'ca 6 'a* a 6 'a* u 

o o o o 



- compare with Equation 26 for the non-uniform flow field - 






2Ca 6 a 
o o 



•)trz^n6* = (or:;,' ( 1+S) - 



'2Ca 6 'a--' 
o o 



Kn-1) , 

— I 2nj 

a* -“o 



( 26 ) 



The transformed Meek Criterion looks very much like equation 26 with some 
sign changes, but the basic form is nearly identical to the non-uniform 
flow field equation. The sign changes can be accounted for by the obser- 
vation that in the non-uniform field equation: 

^0 ~ log(4d/R) 

and as the point radius (sphere or wire to plane equations are similar) 

R (x> ^ the log term becomes negative and the equations merge to nearly 
the same form. 

After seeing how similar the actual physical situations are at 
breakdown though, the similarity should not come as a surprise. In fact, 
what the similarity further suggests is that perhaps there really is no 
absolute uniform field breakdown. In a uniform gap at a potential just 
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below breakdown, statistically, numerous small electron avalanches may 
originate and die out in a short distance as depicted in Figure 10a. Many 
of these mini -avalanches must meet the anode to form minute deformations 
of the field there, any one of which, if large enough, can ultimately 
trigger total gap breakdown similar to the non-uniform point-plane. 

This brings us to an interesting obeservation concerning a highly 
non-uniform field situation where the accumulation of space charge around 
a positive point tends to distort its apparent shape and radius, the result 
of which is to tend to reduce the non - un i form i ty and go toward a more uni- 
form field conf i gurat ion. From these two examples, it would appear that 
the "natural" state at breakdown in a discharge is a moderate state of 
non-uniformity rather than at either extreme. 

6 . Quantitative Comparison of Data and Theory : 

Figures 11 through 17 present some typical data recorded on 
various runs with the point and wi re- to-pl ane electrodes. Conditions 
ranged from high reverse flow velocities (^120 m/s) to high normal flow 
(80-100 m/s with turbulence). No turbulence could be introduced in the 
reverse flow mode, since as illustrated in Figo 1, the cathode is constructed 
of aluminum honeycomb flow straightener which assures smooth laminar reverse 
flow throughout the velocity range tested. 

The data presented were collected over a six week period in 
varying weather conditions and are not corrected for relative humidity and 
ambient temperature effects, both of which affect the properties being 
measured. Correction for these effects might reduce the scatter evident 
in various data sets, especially in the composite summary plots. An 
additional source of statistical scatter due to unequal current flow to 
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various pins or wi res may also be significant, especially in a discharge 
such as ours which is not ballasted to produce uniform current flow from 
each element. As noted in Ref. 3, however, effective turbulence does 
cause more uniform activation of all pins and therefore reduces the 
adverse effects due to no ballast. 

Figure 11 displays corona initiation and breakdown voltages with 
no flow versus gap length for both geometries tested. Theoretical curves 
are also displayed which were computed from the appropriate equations 
developed in Section A . These curves display an intersection of the 
Vb and Vc characteristics at the shorter gap lengths and are very 
similar in appearance to those presented in Reference 8 for single 
point-plane geometries. The intersections define what is termed the 
"existence*' region in Reference 8 for the corona discharge^ To the 
left of this point, a corona discharge will not form, and the gap will 
breakdown (arc) upon reaching the value of Vc . To the right of that 
point a corona will form which enhances the discharge current. A common 
empirical expression which has been found to quite closely fit current 
data for various non-uniform geometries can be written as follows: 

I = KCV(V-Vc) 

The initial constants contain geometric and gas parameters. This equation 
highlights the dependence of I on V and on the difference between the 
applied and corona initiation voltages. The existence regions for the 
present configurations are seen to be for - 

Points-Plane d >- 0.2 cm 



3 ^ 



Wi res-P lane 



d > 2.0 cm. 



Figure 12 presents a comparison of Vb*) from Equation k2 

max 

when 6* 0 . Due to power supply limitations, insufficient points were 

accessible in the wire/plane configuration because of higher breakdown 
voltage requirements to adequately define the wire/plane curve so only the 
point/plane data are presented. This curve is interesting in that it 
shows at longer gaps flow cannot affect breakdown voltage as radically as 
it does in a short gap. The general trends of the Vb* vs. velocity 
data taken under similar conditions of humidity and temperature tend to 
confirm this. However some contradictions are evident and this cannot 
be considered a conclusive result. 

Figures 13, 1^, and 15 are comparisons of data and theory for 
three gap lengths of Vb* vs. velocity . The values of the parameters 
C, m, n, and turbulence index used in the calculated curves are given. It 
has been observed that C = 1.3 , m = 3 , and n= 8 generally result in 
a good first estimate for the basic curve, and then the turbulence index 
(input of turbulence intensity) is varied to match the data. in using 
this technique, it has been observed that long gaps usually require a 
lower turbulence index even though the actual turbulence being generated 
at a given velocity is the same as in shorter gaps. One possible explana- 
tion for this is that, as pointed out by Reference 2 after the turbulencee 
being generated is fully developed, it begins to rapidly decay with 
time and distance, so in a long gap, although turbulence levels entering 
the discharge region are the same as in short gaps, the turbulence 
level leaving the discharge will have decayed more, resulting in a lower 
average turbulence intensity in the longer gap. This would indicate that 
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a weighting of the turbulence index with gap length may be required. 

Probe measurements with screen #9 installed in the test tunnel show that 
at the highest attainable velocities, the turbulence intensity peaks just 
at the beginning of the discharge region (at the pin tips), which would 
support the above argument, and which also may be the key to the success 
of the screen #9 design. 

Figures 16 a, b and 17 a and b are summary plots of Vb and lb 
data vs. velocity for the points and wires. While many gap lengths and 
turbulence levels are included under varying ambient conditions, the plots 
serve to point out the main features and general trends of these parameters 
with flow and turbulence. The lines on these curves are not calculated, 
but are drawn in to highlight the basic features of the various data 
collected and to show the similarity with the theoretical curves dis- 
played elsewhere. 

The next section qualitatively examines the wire/plane in more 
detail. 
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B. Wi re- to-P lane Geometry 

1 . Qualitative Discussion of Differences from Point-Plane 

One of the significant changes in going from the point-plane 
geometry to the wire-plane is a major modification in the shape of the 
electric field strength distribution across a gap of given dimension and 
a fixed applied voltage. While the assumed distribution (proportional to 
X is still “hyperbolic'* as illustrated in Figure 3, at the electrode 
surfaces the point field is much higher and drops off much more quickly 
than a wi re field of the same radius at the same potential. Since in 
both cases the cathode is grounded, the wire electric field in the drift 
region must be higher than that of the point-plane. The wire radii used 
in the current research started at 0.08 cm, but it was found that to 
operate within power supply capabilities, the wire size had to be reduced. 

As a result the wires were incrementally filed down to a wedge shape of 
approximately 30^ angle and maximum thickness of about 0.05 cm. This 
permitted some current enhancement within available power supply limits. 

Several things are evident in comparing point and wire distribu- 
tions (see Fig. 3 )• First, for a given gas, the ionization field strength 

O' 

E. determines the length of the anode extension 5 , so i f E. is lower 
I i2 

than the point of crossing (point P) of the curves at a given voltage, 

5 will be greater than 5 , but if E. is higher than this point, 

w ^ p I ^ 

then 5 will be greater than 5 . Since as will be discussed later, 

p ^ w 

the current is a function of 5 , this behavior can have important con- 
sequences . 

Another important characteristic difference of the point and wire 
distributions lies in the transition region from the high field to the drift 
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field. The transition is well defined for a sharp point, but becomes 
very slow and indefinite in wires and points of large radii. This 
results in a less well defined anode extension in the flatter distributions, 
which should affect the flow interactions in the loosely bound outer por- 
tions of the space charge (6) region. Changes in the behavior of the 6 
region with flow should show up as changes in Vb* with velocity. Some 
major changes in the wire-plane and point-plane Vb* vs. velocity curves 
have been well documented in the current research. Figure 18 illustrates 
the basic nature of these changes. 

It has been consistently noted that the wire curves, while still 
exhibiting the low velocity depression of Vb* , almost always take an 
early (i.e. , low velocity) upturn toward the asymptote whether the flow is 
turbulent or not (Fig. 8). The convective depression, caused by an initial 
lengthening of 6 , is very small, and it is thought that the rapid rise in 
Vb- is due to a rather large movement and dissipation of loose space 
charge at the outer limits of the anode extension. In this region, due to 
the flatness of the E-field, the ion slip parameter (S = u/k.E) can attain 
a significant value over much of 6 at once, with the result that the outer 
portion of the anode extension actually blows off or becomes too tenuous at 
low velocity to act as an extension, which then causes a large and rapid rise 
of Vb* . The asymptotic limits of Vb* in both the normal and reverse 
flow directions are interpreted again as the requirement to break down the 
entire physical gap, as in the poi nt-to-pl ane development. 

A third characteristic of the wire geometry is the inherently 
higher values of electric field in the drift portion of the gap. This 
Increase in the drift field increases the probability that high ionization 
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and streamer activity in the region near the wire will be able to develop 
and sustain self-propagation through the drift field region to the cathode 
by the streamer mechanism discussed in most texts on the subject (see Refer- 
ences 8, 12, 15, 16) resulting in breakdown. On the other hand, however, in 
a discharge designed for optimum energy transfer to a laser gas (pumping), 
the design of the electrodes will to some extent be dictated by this pumping 
requirement. In practice, a trade-off is required between electrode design 
optimized for total power input to the gas and that design optimized for 
best E/N in the drift region for most efficient pumping. In an electric 
discharge laser cavity, for Instance, an input power of 1000 watts/cm^ 
means very little if only two percent is going into the lasing states due 
to an inefficient E/N in the pumping region. 

One of the attractive features of the wi re- to-pl ane geometry is 
the increased ionization volume available continuously along the wire. The 
point-to-pl ane ionization volume, on the other hand, is confined to what 
would amount to discrete points along the wire. 

The point-plane produces discrete 3"dimensional axially symmetric 
multiplication volumes, at each point, while the wire-plane system is a 2- 
dimensional continuous volume spanning the length of the wires. Most elec- 
tron multiplication occurs within these volumes, and so it seems as i f a 
wire should produce more current. But the problem is complicated by the 
electric field distribution which determines the longitudinal dimension and 
a the multiplication coefficient. As will be seen total current is not 
a function of a and 5 alone, but a function of the product (ao) or 
(a*6*) in the flow situation. 

2. Total Current Considerations 

As pointed out in Part II, Section 3, since the total current at 
the anode surface is due to electrons entering the anode from the 
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multiplication region (the positive ions being repelled), the total discharge 
current can be expressed as a function of the number of electrons entering 
and their drift velocity (the assumption is made that if negative ions are 
present, they give up their electron, soon after drifting into the high 
field region; this assumption is supported by Reference 15 and some 
sample calculations presented in Reference 8). 

A single electron entering the high field region will exponentially 
multiply to produce 



S adx 

X 

n = n e 
o 



electrons at the anode surface with n = 1 . Since the current density is 

o ' 

the current per unit area. 



j = en V = i/A (50) 

-'e e e 

the total current (which must be constant across the discharge) can be 
expressed as 

adx 

i = i^e = j^A (51) 

This allows the total current to be discussed in terms of electrons 
arriving at the anode. That number is strongly dependent upon how many 
electrons enter the multiplication region, whether as free electrons, or 
attached to atoms or molecules as negative ions. The electron input is 
described as i^ in the above equation and is 

i =^(en V +env)AH y.e 

o e e - - ‘ 

o 



ko 



(52) 



where 



= free electron density from drift region 
o 

n_ = negative ion density from drift region 
Vg,v. H particle drift velocities 

A = area through which particles pass into multiplication region 

Y = secondary electron supply rate. 



Without covering many details or methods already discussed in 
Section A, the above considerations lead to the normalized current equation 
below 



i 2Ca 6 6“"(a*6*-l) 

i - = = Y*e ^ ^ 

' o 



(53) 



where 

Y* z \ /\ = normalized secondary electron supply rate* 

From equation 52 it can be seen that the units of Y ^re electrons/second, 
and that it is an electron flow rate into the multiplication region* It is 
the counterpart of the Townsend secondary ionization coefficient in a non- 
flowing discharge. 

Understanding the nature of y > ^ able to analyze 

possible flow effects* Electrons coming into the ionization region can 
arrive in two forms, either as free electrons at the characteristic 
electron drift velocity which is generally much higher than typical 
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flow velocities, or attached to a heavier particle in the form of negative 
ions, etc., whose drift velocities and therefore delivery rates can be 
greatly altered by the flow. 



Now we discuss the third component of current, the positive ion 



contribution. In the determination of total current by counting electrons 
flying past the anode surface, the positive ions play an indirect role by 
supplying secondary electrons primarily through impact with the cathode where 
energetic collions can result not only in neutralization of the ions, but 
also release of additional electrons which then begin to move toward the 
anode. These secondary electrons may proceed as free electrons or attach 
to particles and proceed as negative ions, thus the positive ion contri- 
bution shows up as an enhancement of both y and y . Since both 

e n* 

positive and negative ion drift velocities are of the same order of magnitude 
as the flows being discussed, it can be anticipated that the total current 
would be significantly affected by flow through ion slip - and electron 
trapping discussed in Part I lA. 

Considering the free electron component alone, the total input of 
electrons to the multiplication zone is, from equation 52 



At first it would seem that electron velocity in the drift region should 
not be significantly affected by the flow due to low electron slip and 
negligible trapping due to space charges. However 




(5A) 



V, ~ k E(x) 
e e 



(55) 
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!n nearly all non-uniform geometries, the electric field is given by 
expression of the following form - 




cV 

ToglCONSTJ 



with the result that 



E*(x) = V* 



and at breakdown, 



Eb* = Vb* 

Therefore from equations 5 ^, 55 , and 58 - 

i n V 

i* = ^ = Vb-^’^n * 

o I - n V e 

oo eo eo 

This equation shows that even free electron velocities can be altered 
indirectly by flow modifications which affect the breakdown voltage 
and hence the electric field even in regions where electron trapping 
not a significant factor. Vb* and therefore electron speed can be 
by 30^ or more at moderate flow velocities. 

Now the total current can be written as follows - 



( 56 ) 
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i* = n *Vb*e 
e 



Ca 5 6 (a*<S*-l) 

• o o 



The n * term is a function of the total current, and all associated 
secondary electron supply processes such as cathode emission, photoioniza- 
tion, etc., whose variations with flow are as yet undetermined. 

In a discharge where many of the secondary electrons are liberated 

by positive ion impact with the cathode, however, n * will be a function 

e 

of the ion slip and the current equation may assume the following form - 



i* = Vb^ 



(1+S) + 



XI 



Cd 6 5""(a*5*-l) 
o o 



( 61 ) 



Neglecting y* temporarily. Equation 61 will be used to qualita- 
tively examine the behavior of the breakdown current with flow and turbu- 
lence and to attempt to explain some of the observed differences in the 
point- to-plane and wi re- to-plane Vb* and lb vs. velocity data collec- 
ted during the current research effort. 

Figures 16, 17 and 18 are typical plots of breakdown voltage and 
current vs. velocity for the point and wire configurations and will be 
referenced frequently in the following discussion. 

a. Point-Plane Current Interpretations 

Since breakdown voltage variations with flow were already 
discussed in detail in Part II. A. 5, we will concentrate on some typical 
point-plane current data first. Beginning at the origin (U=0) in Fig, 
17a, it can be seen that, as velocity increases in the normal flow 
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(positive) direction, the current begins to rise (region 1), reaches a low 
maximum (region 2) and then declines slightly (region 3)* A typical 
Vb* curve (Fig. 8) shows a depressed breakdown voltage in this region, 
which from Equation 61 tends to decrease lb but the increasing velocity 
term helps to compensate for the reduction as does the product (a*6*) in 
the exponential, and the current in regions 1 thru 3 is dominated by the 
behavior of 6* (in this velocity region turbulence intensity is not yet 
sufficient to cause significant enhancement of a* ). As velocity increases, 
if an effective turbulence generator is in place, a* becomes dominant due 
to its rapid exponential growth with increasing velocity (Equation AO) , and 
in accordance with Equation 61 an overexponential rise in lb is noted 
(region 4 of Fig. 17a). The figure also shows Vb* leveling at the normal 
flow asymptote which would also cause a leveling of the Ib characteristic 
or at least a leaning toward a (1+S) linear rise in accordance with Equa- 
tion 61. 

b. Reverse Flow Current 

Interpretation of the reverse flow (negative velocity) current 
and voltage characteristics in the point-plane and wire-plane is essentially 
the same, so the following discussion applies to both configurations. Since 
no turbulence was introduced in the reverse flow mode, other processes tend 
to dominate breakdown characteristics in this regime. Equation 45 shows a 
rapid decrease in 6* with reverse flow which causes an equally rapid rise 
in Vb* ; the Vb* increase and (1+S) velocity term decrease might be 
expected to neutralize each other as in the normal flow mode, leaving the 
major influence on Ib the 6* term again. Since 6* is rapidly decreas- 
ing, this line of reasoning would have Ib decreasing at a fairly rapid 
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rate. However, the experimental data consistently show a steep rise 
in lb at relatively low negative velocities (region 5 of Fig. 17a), and 
a leveling or decrease at higher velocities (regions 6 and 7). The high 
velocity reduction in lb is easily understood via Equation 61, since the 
velocity term (1+S) will eventually approach zero at high enough nega- 
tive velocities causing lb to approach zero. Also since 6* is con- 
stantly decreasing with increasing reverse flow, this tends to reduce lb 
more at high velocity. (In one experimental case lb was reduced to 
approximately 1/10 of its no flow value). It seems that the only factor 
going in favor of increasing reverse flow current in the low velocity 
region is the increase in Vb* ; but numerical calculations show that the 
velocity term more than offsets this factor. Careful review of the equa- 
tions brings aid from an unexpected source to explain this strong low 
negative velocity current enhancement through the a* term. 

Recalling the form of a from Equation 6, it can be seen that 
a is a function of E as well as P , and as Vb* increases in the 
reverse flow mode, E also increases, since from Equation 1 E - V . 
Mathmati ca 1 ly this leads to the following enhancement of a due to the 
rise of E for the point-plane - 

a*(AE) = exp{(-BP/Vbo)j;(R log(Ad/R)](l/Vb*-l) } (62) 

Another subtle effect may also come into play here, since, as flow is 
directed into the quasi-steady no-flow anode extension region, the initial 
hydrodynamic effect will be to promote a flattening and radial spreading 
of the dense space charge boundary. In this case the point radius R 
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which will be a complicated function of the 



should be replaced by 
flow and electric fields* Qualitatively, however, it is expected that after 
initial low velocity spreading which enhances a* (via Equation 62) and 
therefore i* , i* will eventually decline due to the continued decreases 
in 6- and the (1+S) velocity term in accordance with Equation 6l and 
will also decline. The end result of these processes is a rapid 
reverse flow current rise and fall as illustrated by Fig. 7a., These two 
reverse flow effects a*(AE) and R peculiar to reverse 

flow, but in the normal mode they must compete with other strong influences 
which may tend to mask their effects except under special circumstances. 

For example, in the normal mode, gas flow and positive ion repulsion from 
the anode act in the same direction on the high density positive space 
charge, so gas convection by itself does not promote radial spreading as 
it does in the reverse flow mode where electrostatic repulsion and convec- 
tion oppose each other. If, however, the normal flow becomes highly tur- 
bulent, the random gas velocity and density fluctuations in the ionization 
region tend to enhance ionization and the transport properties of the gas 
resulting in faster radial diffusion and dynamic spreading of the positive 
charges, and as in the reverse flow case, the effective radius of the 
virtual anode is increased. Reference 3 presented photographs of various 
discharge conditions which visually illustrates this turbulent enhancement, 
the glow region in the turbulent case being 3~5 times the radius of the 
no-flow or laminar flow cases. 

c. Wire-Plane Normal Flow Breakdown Current and Voltage 

Breakdown data for this configuration was limited to a short 
gap length range above the existence point for the corona discharge (^2 cm) 
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due to power supply problems which limited available voltage to approxi- 
mately 33KV. This limitation effectively restricted wire-plane gap length 
to 2,8 cm or less with flow. Since this is less than 1 cm into the 
corona existence region, it is not possible to give a full accounting of 
flow effects. In general, however, non-turbulen t curves were usually 
flat with little enhancement of current, while most curves with turbulence 
screen #9 In place showed slow but steady improvement with velocity. One 
interesting test was performed with a gap length below the corona existence 
point, and the introduction of turbulent flow resulted in the flow of a 
small current before breakdown. This is most likely due to the pressure/ 
density fluctuations in the vicinity of the wires which cause corresponding 
fluctuations in E/P . If the lowest pressure variations result in a high 
enough value of the reduced field (E/P), some ionization will begin to 
occur before breakdown. 

The wire-plane breakdown voltage curves were mentioned in Section II 
A. 6 and the general trends shown in Figs. 8 and 18. One additional obser- 
vation will be made here. It has been noted in the point-plane configura- 
tion in the normal flow mode, that as gap length is shortened toward the 
corona existence point, breakdown voltage will rise at lower flow veloci- 
ties. In the most extreme point-plane case, Vb* had reached the normal 
flow asymptote at 40 m/s with no_ turbulence. In fact the curve looked 
almost like a wire-plane curve. This seems to infer that if gap length is 
in the corona existence region, but close to the existence point, the 
normal flow Vb* characteristic may tend toward the wire-plane shape of 
Fig. 18. That is, if the gap is short relative to the corona existence 
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polat a low velocity rise in Vb will occur. Additional data at various 
gap lengths are necessary to confirm this preliminary observation. If, 
in fact, this is a true characteristic, it is most likely related to the 
variation of the axial breakdown electric field distribution as the gap 
length is reduced. 

The remaining sections review requirements for the optimum 
pumping of a laser gas mixture vs. electrode geometry, a final over- 
view of data and theory developed earlier, and some conclusions and 
recommendations for further research. 
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C. Electrode Design/Optimization 

In electric discharges designed for laser pumping, two basic require- 
ments must be satisfiedo First, the electrodes must efficiently produce 
a sufficient electron population to permit effective energy transfer to 
the laser gas. Secondly, the electrode design must be such that the 
reduced electric field (E/N or E/P) in the pumping region optimizes elec- 
tron energy for maximum energy input to the desired vibrational or 
electronic energy levels. Degradation in either of these critical areas 
will result in reduced laser power output. 

Discharge cavity size limitations or specifications will most likely 
dictate the maximum gap spacing between anode and cathode, so within that 
limitation the combination of the stressed electrode radius and gap length 
must be determined which best satisfies the requirements above. Since the 
E-field distribution is a function of both parameters, they cannot be 
considered separately. In addition, it has been shown that space charge 
effects via the “anode extension” and radial spreading cause significant 
changes in the “apparent” electrode geometry which should be taken into 
account in the final design. The gasdynami c interactions with these space 
charge induced dimensional changes, in effect, provides an extra degree of 
freedom in a fixed geometry cavity through some control of the radial 
spreading and anode extension. This in turn permits, to some extent, ad- 
justment of E/N and electron production via flow and turbulence. 
Additional benefits arise through the enhancement of transport properties, 
improving stability while concurrently increasing electron production. 

Calculation of the reduced fields via the adaptation of the hyper- 
boloidal point-plane, and an assumed wire-plane distribution indicates 

^3 



E/N values In the drift regions of the discharges tested slightly higher 
than the optimum range specified in Reference 1, a condition which can be 



remedied by 
an adequate 



proper adjustment of the geometric parameters, provided that 
power supply is available. 
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111. SUMMARY 



A. Data 

The Vb* data for both configurations are quite similar in the 
reverse flow modes, showing a fast rise to an apparent asymptotic limit 
(Figs. 16 a,b). This is presumed due to collapse of the anode extension 
by the gas flow. 

The normal flow mode shows Vb* initially decreasing for all gaps 
and geometries with an eventual upswing which depends strongly on turbu- 
lence intensity, and appears to also depend on gap length relative to the 
corona existence point, with early (low velocity) rise strongly favored in 
the high turbulence and/or short gap cases (Figs. 8, 16 a, b, 1 8) . 

The breakdown current data in the reverse flow illustrated similar 
behavior in the point and wire-plane configurations, starting with a sub- 
stantial low velocity rise followed by a decline at higher velocities. 

The point-plane normal flow current behavior exhibits several features 
under various conditions. The most promient of these are outlined in 
Fig. 17a. Many smaller scale perturbations are evident in some individual 
sets of data, but the high turbulence, higher current curves seem to 
follow the same basic pattern in the point-plane configuration. The wire- 
plane current curves exhibit much more scatter and trends are not as well 
defined as with the point geometry. This is probably attributable to the 
fact that data were not accessible in a wide range of gap lengths, but 
were restricted to gaps close to corona initiation where Vb and Vc are 
not far apart, and therefore current levels are relatively low. 
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B. Theory 

The theory developed in Part il accounts for the major features of 
the Vb and lb curves described above as can be seen in the various 
figures, but many of the analytical terms can be improved and it is 
believed that further study will lead to analytic expressions for the 
various parameters of the theory which have been mentioned. The clear 
physical basis of the theory lends to easy interpretation of various flow 
related phenomena and should prove useful in visualization and design of 
flow-type systems and possibly in trouble-shooting of existing flow system 
problems . 
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IV, CONCLUSIONS AND RECOMMENDATIONS 



A. Conclusions 

1. The basic anode extension model developed herein appears to 
offer a simple approach to the interpretation and calculation of basic 
gasdynamic effects on high pressure corona discharges. 

2. The flow field interactions discussed are basic, so the approach 
used in modelling the positive point-plane electrode system should be 
easily applied to other gap configurations, polarities, gas mixtures, and 
pressure regimes. 

3. A more rigorous evaluation of various terms in the basic model 
should provide better insight into the dominant forces in various flow 
regimes, and most likely will result in definition of some of the small 
scale features not yet well explained. 

B. Recommendations 

1. Continue data collection on the v/i re-pl ane electrode configura- 
tion when a larger power supply is available (50KV, 50ma should be 

suf f i cien t wi th present hardware). 

2. Develop a good analytical expression for the effective electrode 
radius for insertion into the present point-plane computational model. 

3. Include a*(AE) and radial spreading effects in normal flow 
computational scheme. 

4. Test point-plane in reverse flow mode with turbulence at low 
velocities to verify Equation 45. 

5. Develop comprehensive analysis of flow on y* for various types 
of gas mixtures, i .e, , attaching, non-attaching, etc. 



53 




i 



I 



I 



4> 








“O 








O 






• 


JZ 






(J\ 


•i-l 






=tfc 


OJ 


d> 






o 


U 




c 




c 




(U 


-Q 


<u 




(U 


e 






u 


o 


3 




0 


o 


-D 




to 


>- 


i- 








3 




O) 


C 


h- 




c 


o 






•— 


z: 


U 




4-1 








OJ 


£ 


•— 




u 


3 


o 




<u 


c 


c 




c 




a> 




a> 


S 


-C 




0 


3 


a. 






»— 






(U 


< 


•* 


• 


u 






<T\ 


c 




•— 




(U 


0) 


m 






"O 




C 


3 


o 


1- 


(U 


-O 


c 


a> 


d> 


u 


< 


4-1 


i- 


3 




m 


u 


H- 


to 


2 : 


t/1 




c 






-o 




CD 




c 


Q. 


c 




OJ 


, — 


c 




CO 


<D 


(U 




(U 


(U 


4-1 




T3 


4-1 






0 


t/^ 


CD 




u 








4-1 


1 


fO 




U 




u 




(U 


(U 


4-1 






c 






LU 


fU 










2 






Ol 


0 




• 


1 








CO 


U- 






-M 


NW» 




o; 


c 






u 








3 


*o 






O) 


Q. 














u. 




54 





55 



Figure 1. Electrodes and Turbulence Generating Screen 
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(a ) Normal Flow 




CATHODE 



ANODE 



(b) Reverse Flow 
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(c) Cross -Flow 



Figure 2. VARIOUS ORIENTATIONS OF ELECTRODES 
AND TURBULENCE GRIDS TESTED. 
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Figure 3. TYPICAL POINT -and WIRE -PLANE ELECTRIC FIELDS. 




Figure 4. ELECTRIC FIELD AND VOLTAGE DISTRI- 
BUTION BETWEEN PARALLEL PLATES. 
Subscripts = ( I ) No Space Charge 

(2) Moderate Space Charge 

(3) High Space Charge 
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1 . Applied Field 

2. Space Charge Field 

3. Total Field 



Figure 5. (a) Schematic of Field Lines and Equipotentials. 

(b) Apparent Dimensional Changes Due to Space 
Charge. 

(c) Applied, Space Charge and Total Field 
Distributions. 

X and x' are Critical Distances From 
c c 

Anode Where Electron Avalanching Begins 
(E = E. ) With and Without Space Charge 

Effects. 
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Figure 6. BREAKDOWN VOLTAGE RATIO 
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Figure 7. PLOT OF BREAKDOWN VOLTAGE VERSUS FLOW VELOCITY 
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Figure 9 . PLOT OF BREAKDOWN VOLTAGE VERSUS FLOW VELOCITY. 
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SECTION A-A 




Figure 10. (a) Activity in a Parallel Plate Gap. 

(b) Modified E- Field Due to Space Charge. 
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GAP LENGTH (cm) 

Figure 11. BREAKDOWN AND CORONA INITIATION VOLTAGE Vs. GAP LENGTH 
FOR NO-FLOW AT ATMOSPHERIC PRESSURE, FOR POINT- PLANE 
AND WIRE-PLANE GEOMETRIES. 
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FOR POINT-PLANE GEOMETRY. 
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Figure 13. PLOT OF BREAKDOWN VOLTAGE VERSUS FLOW VELOCITY 
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Figure 14. PLOT OF BREAKDOWN VOLTAGE VERSUS FLOW VELOCITY. 
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Figure 15. PLOT OF BREAKDOWN VOLTAGE VERSUS FLOW VELOCITY 
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Figure 16o. POINTS TO PLANE Vb* Vs. U DATA SUMMARY. 
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Figure 16b. WIRES-PLANE Vb"^ Vs. U DATA SUMMARY. 
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Figure l7o. BREAKDOWN CURRENT Vs. FLOW VELOCITY POINT- PLANE 
SUMMARY PLOT. 
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Figure 17b. BREAKDOWN CURRENT Vs. FLOW VELOCITY WIRE-PLANE (3) 
SUMMARY PLOT 
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Figure 18. TYPICAL POINT and WIRE - PLANE Vb* Vs. U DATA PATTERNS. 
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